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Abstract. On any compact Riemannian manifold {M,g) of dimension n, the L^- 

normalized eigenfunctions {4>\} satisfy H^aIIoo < CX 2 where —A4>\ = X?(t>\. 
The bound is sharp in the class of all (M, g) since it is obtained by zonal spherical 
harmonics on the standard n-sphere S" . But of course, it is not sharp for many 
Riemannian manifolds, e.g. flat tori R"/r. We say that 5", but not R"/r, is 
a Riemannian manifold with maximal eigenfunction growth. The problem which 
motivates this paper is to determine the (M, g) with maximal eigenfunction growth. 
Our main result is that such an (M, g) must have a point x where the set Cx of 
geodesic loops at x has positive measure in S*M. We show that if (M,g) is real 
analytic, this puts topological restrictions on M, e.g. only M = (topologically) 
in dimension 2 can possess a real analytic metric of maximal eigenfunction growth. 
We further show that generic metrics on any M fail to have maximal eigenfunction 
growth. In addition, we construct an example of (M, g) for which Cx has positive 
measure for an open set of x but which does not have maximal eigenfunction growth, 
thus disproving a naive converse to the main result. 



1. Introduction 

The problem which motivates this paper is to characterize compact Riemannian mani- 
folds (A/, g) with maximal eigenfunction growth. Before stating our results, let us describe 
the problem precisely. 

We first recall that the associated Laplace-Beltrami operator A = Ag , has eigenvalues 
{— A^}, where < Ag < < A| < . . . are counted with multiplicity. Let {(p^{x)} be an 
associated orthonormal basis of L^-normalized eigenfunctions. If A^ is in the spectrum of 
—A, let Vx = {(f> A(/) — ~X^(f>} denote the corresponding eigenspace. We then measure 
the eigenfunction growth rate in terms of 

(1) L°^{X,9)= sup m\L^. 

l*ll_f,2=l 

If e\{f) is the projection of / € onto V\, let E\f = J2\ <\^x^{f) be the associated 
partial sum operators. Note that 

(2) E),{x,y) = Mx)My) 

A„<A 
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is the kernel, where are as above. If 



(3) pM^^Yj^^^^^ 

is the principal symbol of \/— A, then by the local Weyl law [Av], |Ho 1| 



(4) Ex{x, x) = > ^ \ci>,{x)\' = (2^)-" / + R{\ x) 

Jp{x.,i)<\ 



with uniform remainder bounds 

|i?(A,x)| < CA"-\ x€M. 
Since the integral in is a continuous function of A and since the spectrum o f th e 



Laplacian is discrete, this immediately gives |(/)y(a;)p < 2CA"' ^ (see Lemma 3.4) 



which in turn yields 

(5) i-(A,5) = 0(A"^) 

on any compact Riemannian manifold. 

The bound (||) cannot be improved in the case of the standard round sphere, (5", can), 
(by the zonal spherical harmonics) and on any rotationally invariant metric on 5*^. It is 
however not attained by metrics of revolution on the 2-torus T^, and the problem arises 
whether there is any metric g on a surface of positive genus for which it is attained. 

More generally, we pose: 

• Problem: Determine the {M,g) for which L°°{\^g) = f7(A~2-). 
Here we are using the notation that means 0(A^^) but not o(A^^). 



Our main result. Theorem 1.1, implies a necessary condition on a compact Riemannian 
manifolds (M, g) with maximal eigenfunction growth: there must exist a point x € M 
for which the set 

(6) C^ = {Ce SIM : 3T : exp^ = x} 

of directions of geodesic loops at x has positive surface measure. Here, exp is the expo- 
nential map, and the measure \n\ of a set fl is the one induced by the metric g^ on T*M. 
For instance, the poles xn,xs of a surface of revolution (5^, g) satisfy \Cx\ ~ 27r. 

Theorem 1.1. Suppose that \Cx\ = 0. Then given e > there exists a neighborhood 
M = A/'(e) of X, and a positive number A — A(e), so that 

(7) sup < eA("-l)/^ XespecV^>A. 

If one has \Cx\ =0 for every x € M then given £ > there exists A{s,p), so that 

\\'I>\\lp{M) ,S(v) \ r- I a" \ l^i \ ^ 2(n+l) 

—7, — ^ — -<e\^^' a e spec V— a > a £, n , v > \. , 



(8) sup <eX'^P>, AespecV-A>A(£,p), p> 

where 



(9) S{p) 



"-ICl - M O < -n < 2(71+1) 



RIEMANNIAN MANIFOLDS WITH MAXIMAL EIGENFUNCTION GROWTH 



3 



(10) 



The LP-bonds improve on some estimates in |So2|, which say that 

<P\\p 



sup 



2 < p < oo. 



The final result is false in the endpoint case p = '^'"^^P , as will be explained in ^ 

Recently, Taylor | Ta | proved bounds lik e ([z] ) under more restrictive hypotheses involv- 
ing the caustics of the wave group, and in |Ta he also discusses applications of these sort 
of bounds to problems involving the convergence of eigenfunction expansions. Also, Sar- 
nak |Sa| has recently obtained limit formulae for the L^-norm of certain cigenfunctions 
on arithmetic hyperbolic quotients. 

Our eigenfunction bounds are based on the following estimates for the local Weyl 
remainder. 

Theorem 1.2. As in let R(X,x) denote the remainder for the local Weyl law at x. 
Then 

(11) R{\x)^o[X'-^) if |/:,| =0. 

Additionally, if\Cx\ ~ then, givene > 0, there is a neighborhood M of x and a A =< oo, 
both depending on e so that 

(12) \R{X,y)\<eX^-\ yeN, X>A. 

As we were informed by V. Guillemin, A. Laptev and D. Robert after completing this 



paper, these estimates overlap previous results of Y. Safarov in the article [sai] (see also 
pV[ , §1.8 and |[v2|). Because our proof is somewhat different from that of Safarov and 
our exposition contains complete details, we have decided to retain our original proof. 



Theorem (1.2) is a kind of local analogue of the classical Duistermaat-Guillemin [DG| 
(see also Ivrii[|[vl|) result that N{X) = (27r)^"Vol(M)A" + o(A"-i) if the set of initial 
directions (x,^) of closed geodesies has measure zero in S*M. Here, N{X) denotes the 
number of eigenvalues < A (counted with multiplicity). Their assumption is the natural 
one for this theorem since only closed geodesies contribute singularities to the trace of 
the wave group U{t), t 0; however, for pointwise bounds at x one needs to assume zero 
measure of closed loops through x. 

We should point out that ( [Til ) by itself is not strong enough to give the above sup- 
norm estimates for cigenfunctions. They also require the more delicate second statement, 
which says that x limsup^;^^;^ A"~^i?(A, x) is continuous at points where \£x\ = 0. 
We now give some concrete applications and examples of these general results. 

Our first application is to the real analytic setting, the result has strong implications 
on the topology of M. 

Theorem 1.3. Suppose that {M,g) is real analytic and that L°°{X,g) = ri(A("^^^/^). 
Then {M,g) is a -manifold, i.e. a pointed Riemannian manifold {M,m,g) such that 
all geodesies issuing from the point m return to m at time £. In particular, if dim M = 
2, then M is topologically a 2-sphere . 



For the definition and properties of -manifolds, we refer to |Besse| (Chapter 7). 
By a theorem due to Berard-Bergery (see |BB, Bessej, Theorem 7. 37), manifolds 
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M satisfy 7ri(Af) is finite and H*{M, Q) is a truncated polynomial ring in one generator. 
This of course implies M = S"^ (topologically) when n = 2. We remark that the loops are 
not assumed to close up smoothly. An interesting example to keep in mind here is the 

2 2 2 

tri-axial ellipsoid -Eai,a2,a3 = {{xi,X2,xz) € : + + = 1}, with ai < 02 < 03. 

This is a metric on S"^ with four umbilic points. At each umbilic point, all the geodesies 
leaving the point return to it at the same time. Only two such geodesies close up smoothly, 
namely the middle closed geodesic {x2 = 0}, traversed in the two possible directions. For 
more on the geometry of the ellipsoid, we refer to @, p^, [CVV[ . 

Our second application is to generic metrics: 

Theorem 1.4. L°^{X,g) = o(A("^i)/2) for a generic Riemannian metric on any mani- 
fold. 

The proof is just to show that the condition \Cx\ = 0, Vx G M holds for a residual 
set of metrics with respect to the Whitney C°° topology. It appears that this is a new 
geometric result. 

To our knowledge, the most general class of metrics for which an improvement of the 
general error bound 0(A~2~) in the local Weyl law has previously been proved is that of 
manifolds without conjugate points satisfying exponential bounds on the geodesic flow 
In that case, the error term has been improved to 0{ ^ ^^^^ — ). The result is not 



stated explicitly in [ Be ] , but follows from its estimates on the remainder term in the Weyl 
law. Estimates of the error term in the Weyl law also appear in Q , but it appears that 
the methods of that paper require integration in x and hence do not imply bounds; 
in particular, their assumptions are on closed geodesies rather than loops. 

It is natural to ask whether the converse of these results is also true. Regarding the 



first remainder estimate in Theorem 1.2, we shall prove the conditional converse: 

Theorem 1.5. If the loopsets of{M,g) are clean, or even almost-clean, then the follow- 
ing implication holds: 

R{X,x)^o{r-^) =^ \Cx\ ^0. 

Clean and almost-clean loopsets are generalizations of the same notions for fixed point 
sets and will be defined in §0. Almost-clean (M, g) include manifolds containing 'round 
annuli', i.e. annuli isometric to equatorial annuli of the round sphere. One can obtain 
such manifolds by deforming the round metric on S" in small balls, say by adding handles. 



The geodesic flow in such cases is not clean in the sense of Duistermaat-Guillemin [DG| 



but is 'usually' clean outside a hypersurface (the boundaries of the balls). The proof of 



Theorem 1.5 is similar to the proof of the converse to the Duistermaat-Guillemin theorem 
in jz| and the assumption of almost-cleanness is taken from there. 

However, the naive converse to our sup-norm result is simply false. 

Theorem 1.6. There exist C°° -Riemannian tori of revolution (T'^,g) such that: 

3x:|£,|>0, i?(A,x) =17(A^), but ^^^^^ = o{X^). 



Thus, remainder blow-up is not sufficient for maximal eigenfunction growth. We note 
that it is very likely that the ellipsoid is also an example of such behaviour: we have 
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i?(A, x) = ri(A"2 ^ ) at each umbilic point, but it is conjectured that — o(A 2 ) (J. 

Toth). In the final section, we shall discuss some interesting open problems on maximal 
eigenfunction growth. 

This paper is organized as follows. In ^ , we review some basic results on the wave 
group and geodesic flow, a nd o n se vera l Tauberian theorems. This back ground is sufficient 
for the pr oofs of Theorems LI and L2 in §|[ The proof of Theorem 1^ is in par t a va riant 
of Ivrii's |Ivl] (see also [ Ho IV| proof of the Duistermaat-Guillemin theorem DG | that 
the remainder term in Weyl's law is o(A^~) if the set of closed geodesies has measure 
zero). In particular, we shall need to analyze the small-time behavior of the restriction to 
the diagonal in Af of the kernel of CU{t)B, if C and B are zero-order pseudodiferential 
operators. The case where C = Id is the identity operator was analyzed by Ivrii|Ivl| 
[Ho IV I; however, since we cannot take traces in our proof, we need to study this slightly 
more technical case. Fortunately, we are able to reduce the analysis to the special case 
discussed in §|]. Additionally, Theorem 1.2 uses some continuity properties of \Cx\- The 



proof of Theorem 1.1 for the L 



norms is then a simple matter, and the other cases 
We then apply the results to some 

" In§| 

to real analytic manifolds and generic 



follow by interpolating the bounds in [Sol 

concrete classes of metrics. As preparation, w e di scuss the geometry of loops in 
and we give the applications of Theorem 

Riemannian manifolds, respectively. In §^ we study converses to Theorems LI and |l.2| , 
and construct an example where R{X,x) — r2(A"^) but L°°{X,g) — o(A"t~). Finally, in 
we consider related open problems. 

Regarding notation, C will denote a constant which is not necessarily the same at 
each occurrence. Also, integrals like the one in (^) should be interpreted as C°° densities 
which are homogeneous in A and defined so that the scalar product with ip g C°°{M) is 
the integral of the product by -0 with respect to the symplectic measure dxd^ on T*M. 

It is a pleasure to thank R. Hardt, W. Minicozzi, B. Shiffman for very helpful conver- 
sations on analytic sets, and for giving us copies of the unpublished articles ^ |s|. We 
are particularly grateful to V. Guillemin, A. Laptev and D. Robert for the reference to 
Safarov's work [Sal[. We also thank J. Toth for bringing the ellipsoid example to our 
attention. 



2. Wave group and geodesic flow 



In this section, we collect the relevant background on the wave group and geodesic 
flow for the proofs of Theorems L2 and |l.l[ Essentially, they are based on the short- 
time asymptotics of the wave equation, long time wave front relations and Tauberian 
theorems. We shall need substantially more geometry of loopsets to prove Theorems 1.5 
and 1.6, but postpone the discussion until §4 



Throughout this paper, (M, g) denotes a compact Riemannian manifold of dimension 
n > 2, T*M denotes its cotangent bundle and S* AI denotes its unit sphere bundle 
with respect to g. We denote by exptHp the geodesic flow of g, defined as the flow 
of the Hamiltonian vector field Hp of p(a;,^) = v^ff^H^zOOCfc? the principal symbol 
for \/—A. By definition, exptHp is homogeneous, i.e. commutes with the natural R^"- 
action, r x (a;,^) — {x,r£^) on T*Af\0. We also define the exponential map at x by 
exp^ ^ = TroexptHp{x, ^). These definitions are standard in microlocal analysis but differ 
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from the usual geometer's definitions, which takes = ^ 9^'' {^)Ci^j ^-s the Hamihonian 
generating the geodesic flow. The geometer's geodesic flow is not homogeneous. 

2.1. Wave group and spectral projections. Let dE\ be the spectral measure as- 
sociated with the partial summation operators E\ for \/—A. Then we shall prove our 
estimates for the local Weyl law remainder term i?(A, x) by studying the Fourier trans- 
form 

U{t) = I e-'*^dEx = e"^*^^. 



This of course is the wave group, i.e. the unitary group generated by ^/—A. It solves the 
Cauchy problem 



■idt 

It is well known that the kernel U(t, x, y) is a Fourier integral operator in the class 
X M X M,C) where C is the La grangean 

{{t, X, y, T,£„T]) : T+ p{x, 0=0 and (x, £,) = exp{tHp{y, 77))}. 

The last condition means that the orbit of Hp which is at {y, 77) at time is at (x, £,) at 
time t. In particular, x and y are joined by a geodesic of length \t\. 

We shall mainly be concerned with the restriction of the wave kernel to the diagonal 
in M X M. The restriction is defined by 

U(t,x,x) ^ A*U{t,x,y), where A: Rx M ^Rx M X M is A{t,x) ^ {t,x,x). 

Although the notation conflicts with that for the Laplacian, it is standard and should 
cause no confusion. By general wave front set considerations, one always has 

WF{U{t,x,x)) c Ca, 

where 

(13) Ca = {it,T,x,^- T]) : T = - |^|, expiiJp(x,7y) = (x,^)}. 

Indeed, the canonical relation underlying A* is given by (| pG[ |, (1.20)), 

WF'{A) = {{t,T,x,{^ + r]));{t,T,x,^,x,rj))} C T*((M x M) x [R x M x M)). 
The intersection defining Ca is essentially that between C and the second component 
V := {{t,T,x,£„x,ri)} cT*{Rx M x M) 



of WF'{A). If C n X> is a clean intersection (cf. ^G|, |Ho IV[ |, §??), Ca is a Lagrangian 
submanifold of T*{M. x M)\0, and one has the much stronger statement that 

U{t,x,x) e /°(M X M,Ca), 

where /°(]R x M,Ca) is the class of Fourier integral operators of order zero associated to 
the canonical relation Ca- However, we shall not need this hypothesis in the proofs of 
Theorems |l.2Nl.l|. 



We note {x,ri) is the initial (co-)tangent vector and {x,£^) is the final (co-)tangent 
vector of a geodesic loop of length t at x; loops need not be smoothly closed. If there are 
no geodesic loops of length |to| through xq e M, it follows that {t,x) U{t,x,x) is free 
of singularities at {to,XQ). Since the length of a geodesic loop is at least the injectivity 
radius S of {M,g), it follows that U(t,x,x) is smooth for \t\ < S. 
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2.2. Small-time asymptotics for microlocal wave operators. In the proof of Theo- 
rem 1.2, we shah need to calculate asymptotics for the small-time behavior of the diagonal 



restriction [U{t)Q](x, x) {{x,x) : x M} of the one-parameter family of operators 

U{t)Q ; C°°(M) C°°{M), 

where Q{t,x, D^) is a polyhomogenous zero-order pseudodifferential operator. In th e 
case where Q i s in dependent of t, this calculation was carried out by Hormander Ho IV |, 
following Ivri! |lvl| , and, as we shall see, the slightly more general case mentioned above 
follows in a straightforward way from this special case. 

To set up the notation, write 

(14) Q{t, X, D,) = B{x, D,) + tR{t, x, D,), 

where of course B{x, D^) — Q(0, x, D^)- We then shall let h and denote the principal 
and subprincipal symbols of B, while r shall denote the principal symbol of i?(0, x, Dx). 
We recall that p = y^^^pH^^Wj^ denotes the principal symbol of \/— A, and that = 
for this operator. 

The calculation we require then is the following 

Proposition 2.1. For any compact Riemannian manifold (M,g), the restriction K of 
the kernel ofU{t)Q to R x {(x,x) : x G M} is conormal with respect to {0} x {{x^x) : 
X £ Af } in a neighborhood of this submanifold. Moreover, there is a S > so that when 
\t\<S 

(15) Kit,x)= rM^e-^*dA, 



dX 



where A G S"', A{x,0) = and 



(16) A(x,A)-(2^)-" / (6 + 6^)de + (2^)-"|- / i{b,p}d^ 

Jp{x.,^)<\ C*-^ Jp{x,i)<X 



+ (2^)-"^ / ir(0, X, e 

h(x.,o<\ 



dX 



li Q — B and thus is independent of t then (16) is just Proposition 29.1.2 in plo IV |. 
Since the kernel Kit, y) in this case is conormal, it is easy to see that this special case 
also yields the time-dependent case. The constant S is the injectivity radius of (M, g) 



Indeed if 

(17) Q(i, X, D.,) = tR{t, X, Dx)^tY, J\^^^^ 



then we can compute the contributions of the various terms of the expansion to (p^). 
We first note that, by the aforementioned time-independent case, the j = term will 
contribute a term of the form 



Ko{t,x)^ I t^(x,\)e-'^'dX= / -i^e-'^'dX 
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with Ao e S'" satisfying 

Ao(a:,A)-(2^)-" / {r + r^)d^+{2i,)--^ [ l{r,pKe5"-2, 
Jp<\ OA Jp^x 

if r and r** are the principal and subprincipal symbols of R{0, x, D^)- 
Note that 

and so Kq exactly contributes the last term in the right side of ( [l6| ) . Since this argument 
also shows that the remaining piece 

only contributes a term dA/d\ to ( |l5| ) with A G S*"^^, the proof is complete. □ 



To prove Theorem 1.2 we shall have to write 

V{t) = CU{t)B 

as an expansion like (]l5|)-(p^) if C(x, D) and B{x, D) are zero-order polyhomogeneous 



pseudodifferential operators on M. In Ivrii's |[vl| proof of the Duistcrmaat-Guillcmin 



|DG] theorem (see also [Ho IV]) it sufficed just to compute the asymptotics for UB since 
the ultimate calculation involved a trace over M, rather than a pointwise estimate like 
the one we require. 

F ortu nately, though, this slightly more complicated situation follows from Proposi- 
tion To sec this, we note that V{t) solves the Cauchy problem 

\{\§l + P)V{t)^[P,C]U{t)B 
\v{Q)^CB. 

Consequently, by DuHamel's formula 

(18) CU{t)B ^U{t)CB + i f U{t){U{-s)[P,C]U{s)B)ds. 

Jo 

If we change our notation a bit and let aprin{C B) and (7sub{CB) denote the principal 
and subprincipal symbols of CB then Proposition tells us that for small \t\ we can 
write the restriction to the diagonal of the kernel of U{t)CB as 

9Ao, _,xt 



dX 



ix,X)e-'^* dX, 



where 



Ao{x, A) - (27r)-" / {ap„n{CB) + a,ub{CB))d^ 
Jp<\ 



+ ( 



2^^)"" / h{^prrn{CB),p}di G 5"-^ 
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To perform the same calculation for the last term in (|18|), we note that by Egorov's 
theorem 



Q{t,x,D^)=i [ U{-s)[P,C]U{s)Bds 
Jo 



is as in Proposition 2.1 with Q{0, x, Dx) = and dtQ{0, x, D^) = i[P,C]B. Thus, for 
small |t| the restriction to the diagonal of the kernel of the last term in ( p^ can be written 
as 

^^^-{x,X)e-'^' dX, 



Ai{x,X) ~ (27r)-"^ / aj,„n{[P,C]B)dC G S" 



dX 

where 

« r 

^n— 2 

/p<A 

Let b and c denote the principal symbols of B and C. Note then that 

^{ap„n{CB),p} ^ ap„ni[P,C]B) = ^{c{b,p} ~ b{c,p}) . 

Thus, we can combine the main term for Ai with the last term for and obtain the 
following result which will play a key role in the proof of (|l^) . 

Proposition 2.2. The restriction K of the kernel of CU{t)B to K. x {(x, x) : a; G A/} is 
conormal with respect to {0} x {{x,x) : x G M} in a neighborhood of this submanifold. 
Moreover, there is a S > so that when \t\ < 5 

(19) Kit,x)^ r^-^^e-^UX, 



dX 



where Acs G 5"', ^cb(0. A) = ana 



(20) Acb{x, X) - (27r)-" / (c5 + asub{CB)) d^ 



2.3. Tauberian Lemmas. In this section we shall collect the Tauberian lemmas that 
we need. 



The first one is a special case of Lemma 17.5.6 in [Ho III]. It is a slight variant of 
the one used by Ivrii[[vl in his proof of the Duistermaat-Guillcmin theorem. It requires 
a monotonicity assumption that will only be fulfilled for the "diagonal" terms of our 
approximation to Ex{x,x). 

Lemma 2.3. Suppose that ^ is a non- decreasing temperate function satisfying /i(0) — 
and that v is a function of locally bounded variation such that v{Q) = 0. Suppose also 
that TO > 1 and that (p G 5(R) is a fixed positive function satisfying J (j){X)dX = 1 and 
(t>{t) = 0, t <^ [-1, 1]. If (f)cr{X) = a^^4>{X/a), < (J < (To; assume that for XeM. 

(21) \dv{X)\ < (Ao(l + |A|)" + A,il + |A|)™-i) dA, 
and that 

(22) \{{dfi - di^) * (/),)(A)| < B{1 + |A|)-2. 
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Then 

(23) |m(A) - v{\)\ < C^{Aoa{l + |A|)™ + + |A|)™-i + S), 

where Cm is a uniform constant depending only on ctq and our m > 1 . 

The other lemma that we require will allow us to handle the "off-diagonal" terms 
in our approximation to E\{x,x) where the above monotonicity assumption will not be 



valid. This lemma is exactly the one employed by Hormander |Ho 1| and Avakumovic 



|Av] and Levitan [f^ ] in their pro ofs of the Q(A" ^) remainder estimates for the Weyl 



formula. A proof can be found in Ho 1| (or |So3|, p. 128). 



Lemma 2.4. Let g{X) be a piecewise continuous function o/R. Assume that for some 
m > 1 there is a constants Aq and Ai so that 

(24) |g(A + s)-g(A)| <Ao(l + |A|)™+Ai(l + |A|)™-\ < s < 1. 
Suppose further that for some fixed S > 

g{t) = 0, \t\ < S. 
Then there is a constant Cm,s, depending only on m and d, so that 

(25) |g(A)| < Cm,s{Ao{l + \X\r + A,{1 + \X\y 



We also quote a simple asymptotic result of |zj (Lemma 3.6), which uses this Tauberian 
Lemma and which will be used in ^ In the following, SV for any homogeneous V C T*M 
denotes the s e t of un it vectors in V. For background on Maslov indices, excesses etc. we 



refer to ^G|] |Ho IV | 



Lemma 2.5. Let A G L'^{M x Af, A) denote a 0th order Fourier integral operator associ- 
ated to a homogeneous canonical relation A C T* M x T* M , and assume that An Aj^.m 
is a clean intersection. Then 



E 

A„<A 



(A</)„ = z™A^ / a^d/io + 0(A^~i) 

JSAa 



where d/ig is a canonical density on the intersection An /S.t'M) is the symbol of A, 
m is the Maslov index, and e ~ dim A D At* m ■ 



3. Proofs of Theorems [li] and [L2 



We first concentrate on our Weyl formula remainder estimates in Theorem 1.2 and 
postpone proof of our eigenfunction bounds in Theorem 1.1 until the end of this section. 



3.1. Local Weyl law remainder: Proof of Theorem 1.2. Following Ho IV ], we 
begin by introducing the loop-length function on r*Af \0 given by 



(26) 



L*(x,^) = inf{i > : exp^ = x}, 



where L* is defined to be -l-oo if no such t exists. It is homogeneous of degree zero, so it 
is natural to consider the restriction of L* to S*M = {{x,S,) : ~ ^i- ^ 

fact for us is that L* is a lower semicontinuous function, or equivalently that the function 
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1/L*{x,^), which is defined to be zero when L*{x,^) ~ +00, is an upper semicontinuous 
function on S*M. For fixed a; G M we define the set of loop directions at x by: 

(27) C, = {i&SlM:l/L*{x,i)^Q}. 
The compfimentary set 

S*M^\C^ = e SIM : l/L*{x,0 = 0} 
is the set of all unit vectors for which there is no geodesic loop with initial tangent vector 

Let us first show that if xo € Ad and \Cxo \ = 0, then given £ > 0, we can find a ball B 
centered at xq and a A < 00 so that 

(28) \R{\x)\ < eA"-\ x£B, A > A. 

To do this, choose a coordinate patch fi, with coordinates k,{x) around xq which we 
identify with an open subset of M". Let then f{x,^) denote the image of 1/L*{x,^) in 
the induced coordinates for: {{x,^) G S*M : x G fl}. We then need the following 

Lemma 3.1. Let f be a nonnegative upper semicontinuous function on O x S"'~^, where 
O C M" is open. Fix xq € O and suppose that € S"~^ : f{xo,^) 7^ 0} has measure 
zero. Let e > and 6 > be given. Then there is a neighborhood J\f of xq and an open 
set flh satisfying 

\f{x,o\<s, {x,o&MxS"-\nt, 

\nb\ <e. 

Furthermore, there is a 6(^) G C°° supported in flf,, satisfying <b < 1, and having the 
property that if B{^) = 1 — b{^) then |/(a;,^)| <6onAfx supp B. 

Proof. By assumption, 

Es = {^€S^-' : f{xo,0>S} 

satisfies \Es\ =0. Let 

EsU) = e 5"-' : f{x,0 > S some x € B{xo, 1/j)}, 

if B{xo, 1/j)} is the closed ball of radius 1/j about xq. Then clearly, Es{j + 1) C Es{j). 
Also, if ^ G r\j>iEs{j), then for every j one can find Xj G B{x, 1/j) so that f{xj,^) > 5, 
which by the upper semicontinuity of / in turn implies 

6 < limsup/(a;j,0 < /(a;o,0- 

i— »oo 

Thus, 

nj>iEs{j) c Eg. 

Consequently, if j is large, |i?5(j)| < e/2. Fix such a j = jo and choose an open set flf, 
satisfying Es{jo) C fib and \flb\ < e. Then clearly |/(x,^)| < 5 if (x,^) G B{xo, 1/jo) x 

For the last part, note that the argument we have just given will show that the sets 
Es{j) are closed because of the upper semicontinuity of /. Thus, if Es{jo) and f2b are 
as above, we need only apply the C°° Urysohn lemma to find a smooth function b{^) 
supported in Qb with range [0, 1] and satisfying b{^) = 1, ^ G fib, which then will clearly 
have the required properties. □ 
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It follows from this Lemma that if \Cxo \ — 0, then given < £ < oo there cannot be 
too many geodesic loops of length < £ through points x near xq . 



As in Ivrii's proof [|v2| (see also |Ho IV]) of the Duistermaat-Guillemin theorem, we 
now construct a microlocal partition of unity which is adapted to the loop geometry of 
{M,g). It will suffice to construct two pseudodifferential operators b{x, D), B{x, D) — 
I — b{x, D) with the property that B is microsupported in the set where L* {x, ^) >> T. 
We then study E\{x^x) by writing it in the form: 

(29) Ex{x,x) = [{B + b)Ex{B + b)]{x,x). 

The expert reader will observe that this expression is more complicated than the one 

Ex^ B*BEx + b*bEx 



used in |Ivl] |Ho IV |. However, the simpler form is not useful in our problem because 
we are not taking the trace of E\ and cannot move the adjoints in front of E\ to make 
the individual terms into positive operators. Rather, we shall have to expand ( |29| ) and 
study each time of term separately in Lemmas 3.2- 3.3. Moreover we shall have to define 
b carefully to obtain the continuity properties of -R(A, x). We now turn to the details of 
the construction. 



If T > 1 is a large parameter, then by Lemma 3.1 we can find a function b G C°°(S'" ^) 
with range [0, 1] so that 



and 

\f{x,C)\<l/T, onTVxsuppB, 
where N C k,{Q) is a neighborhood of k{xq) and 

BiO - 1 - 6(0- 

Choose a real- valued function ■;/; g C°° which vanishes outside of N and equals one in 
a small ball centered at h{xq). Using these functions we get zero-order pseudo-differential 
operators on R" by setting 



6(x, D)f{x) = i>{x) (27r)-" j j e'^^~y'><b{S,/\i\)i^{y)f{y) d^dy, 

and 

B{x,D)f{x)=^P{x)i2n)-- f f e^^'^-yyi B{^/\^my)f{y) d^dy. 



Note that both variables of the kernels of these operators have compact support in JV. If 
we let b{x, D) and i?(x, D) in '^^{M) be the puUbacks of 6(x, D) and B{x, D), respectively, 
then 

(30) b{x,D) + B{x,D) = ^P'^{x), 

where we are abusing the notation a bit by letting ■0 denote multiplication by the puUback 
of our R" cutoff function. 

Note that by construction 

(31) B{x, <^) = if (x, x, t; r, ??) e C for some {t, r, 77) with < |t| < T, 
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with C being the Lagrangean of U{t, x, y). Consequently the kernels of U{t)B* and BU{t) 
are smooth at the diagonal when < \t\ < T. 

Next, let bo, ¥ and i?o, denote the principal and subprincipal symbols of b{x,D) 
and B{x,D), respectively. Then 

(32) [ Mx,0\'d(+ [ \Bo{x,0\'d^<C, 

where C depends on our local coordinate transformation above but is independent of the 
large parameter T. Also, one has bounds 

(33) / i\bo\ + \Bo\mbo,p}\ + \{Bo,p}M 
Jp{x.,i)<i 

+ / {\as^ABB*)\+asubibb*)\ + 

^ sup 

{Bb*)\)did£_ < Ct, 

Jp{x.,i)<i 

where this constant does depend on T. 

We shall use these facts and the results of the last two sections to prove that for A > 1 

(34) \ij*ix){Exix,x)^i27Ty [ d^)\ <CT-^X"-^ +CtX"'\ 

where Ct depends on T but C does not. This of course yields (|2^). 

In view of (|3^) this bound is a consequence of the following two lemmas. 

Lemma 3.2. There are constants C and Ct as in (^ ) so that if B — B{x, D) and B* 
is its adjoint then for A > 1 

(35) \BExB*ix,x)-i27T)-^ f i\Bo\^ + a,^,iBB*))d^ 

+ {2n)-"dx [ ^{Bo{B*,p}-B*{Bo,p})d^\ 

and if b = b{x, D) and b* is its adjoint then for A > 1 

(36) \bExb*{x,x)-{2n)-" f {\bo\'' + asub{bb*))d^ 

Jpix.i)<\ 

+ {2n)-"dx [ Ubo{bo,p}-b*o{bo,p})d^\ 

Jp{x.,i)<\ 
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Lemma 3.3. There are constants C and C't as in ( |3^ ) so that for A > 1 
(37) 

\BExb*{x,x) + bExB*{x,x)-{27T)-'' f {Bobl + boB*+a,ub{Bb*) + asub{bB*))d^ 

Jpix.,i)<\ 



+ {27r)-^dx / ^{Bo{b*,p}-b*Bo,p + bo{B*,p}-B*{bo,p})d^\ 

One obtains from these three inequahties after noting that ( ^o|) imphes that 

{Bo + bo){B;+b*)=^/j\ 

and 

(Jsub{BB*) + asubibb*) + <7sub{Bb*) + asub{b*B) = asub{i^^) = 0. 

Proof of Lemma ^ To prove (|3|) we note that by (|l|) K{t,x) = BU{t)B*{x,x) is 
smooth for < |t| < T. Therefore if we let v{X) = AsB-'ix, X) be as in (pO[), then since 
( p^ holds for small values of t if /i(A) = BE\B*{x, x) then we have 

|((d^-di.)*<^,)(A)| <i?Ar(l + |A|)-^ 

for constants Bn independent of if < 1/{2T) is fixed. By (||) and (||), (|l|) must 
hold with Aq = C and Ai = Ct where C and Ct are as in (^^. Since 

A„<A 



is non-decreasing Lemma 2.3 yields (|35|). 

The proof of (|3^) follows from the same lines except that here one takes a to be 
twice the injectivity radius of M and uses the fact that Jp^i l^opf^C ^ C/T"^ to get the 
0(T-2A"-1) + 0(A"-2) bounds. □ 



Proof of Lemma 3.3| In this case since 

A BExb* {x,x)+ bExB* (x, x) 



is no longer non-decreasing the first Tauberian lemma, Lemma 2.3, does not apply. In- 
stead we must use the second one, Lemma 

To this end let g(A) ~ gx{X) be given by 

g{\) = {BExb* {x, x) + bExB* {x, x)) - {A*b,{x, A) + AtB- {x, A)) , 



where Asb* and Ai,b* sue given by (20). By Proposition |2.2| , the Fourier transform of dg 
vanishes for \t\ < S, and so supp tg{t) C {0} U {t : \t\ > 6}. But then if {0} G supp g, g 
would include a sum of derivatives of delta functions, which is impossible since g is clearly 
a picccwisc continuous function, and so g(t) = for \t\ < 6. Therefore, by Lemma [2.4| we 
would have ( |37| ) if we could show that there are constants C and Ct as in ( ^^ so that 

(38) |5(A + ,s) - .g(A)| < CT-\l + |A|)"-i + Ct{1 + IMY'-^ 0<s<l. 



Let us first see that this bound holds when g is replaced by Asb* or Ai,b* ■ Here, we 
shall just use the fact that, modulo terms of order n — 2, both Asb-' and AhB-' are the sum 
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of a function which is homogeneous of degree n and one of degree n—1. Consequently, the 
bound will just follow from the size of the coefficients of these two homogeneous terms. 
For Asb* we note that, by (|20|), 

ABb' (x. A) = c(x)A" + d(x)A"-i + 0(A"-'), 

where 



c(x) = (27r)-" / BobodC, 



and 



d{x)^(2iT)-^ f asub{Bb*)d^~n(27r)-^ f ^{Bo{bo,p} - bo{Bo,p})d^. 



By (|3^) and Schwarz's inequality, \c{x)\ < C/T, while (|3^) yields < Ct and so 

( ^ ) must hold when g is replaced by Asb* ■ 

Clearly the same argument applies to Abs-'- Therefore, since 
\BEx+sb* (x, x) - BExb* (x, x) \ = \bEx+sB* (x, x) ~ bE^B* (x, x)| 

= I J2 BM^MAxjl, 

X<X^<\+s 

the proof would be complete if we could show that 

I J2 BMx)bM^\ < CT-\l + \X\r-' + Ct{1 + \X\r-^ < s < 1. 

A<A„<A+s 

To prove this we apply the Cauchy-Schwarz inequality to get 
I BMx)bM^\ 

X<X^<X+s 

X<X„<X+1 A<A„<A+1 

= \BEx+iB*ix,x) - BExB*ix,x)\'/^ \bEx+ib* ix, x) - bE^b* ix,x)\'/^ . 

By ( ^5|) and the argument just given for Asb* 

\BEx+iB*{x,x) - BExB*{x,x)\ < C(l + jAI)"-! + Ct(1 + {XIT'^ 

since the coefficient of the term in the left side of (^) which is homogeneous of degree 
n is bounded by a uniform constant C, while that of the {n — 1) degree term is bounded 
by a constant Ct- 

To finish we need to see therefore that 

\bEx+ib*ix,x) - bExb*{x,x)\ < CT-^{1 + |A|)"-i + Ct(1 + |A|)"-2. 

but this follows from (Bfl) since the coefficient of the term which is homogeneous of degree 
n is bounded by CT^. 



This completes the proof of Theorem 1.2. 
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3.2. Applications to eigenfunctions: Proof of Theorem 1.1. It is now a simple 
matter to complete the proof of Theorem We begin with the following simple fact. 

Lemma 3.4. Fix x G M. Then if X G spec \/—A 

(39) sup ^ = Vi?(A,x-)-i?(A-0,x), 

if R{X — 0, x) denotes the left-hand limit of R at A. 

Proof. Let 

K\{x,y) = ^ (p^{x)(l)^{y) 

A„=A 

denote the kernel of the projection onto Vx- Then if a; G M is fixed, y Kx{x, y) is in 
V\. Consequently, by orthogonality, 

^ \Mx)\'f^^{f \Kx{x,y)\'dyf^ ^ sup | / Kx{x,y)<l>{y) dy \ 



6VA,it'*l|2 = l 

sup \(l>{x)\. 

i>evx, ||0||2=i 



Note that A E\{x,x) is an increasing right continuous function, and so by (^) 

^ =i?(A,a;)-i?(A-0,x), 

A„=A 

which along with the last identity gives ( ^9|) . □ 



The first bound in the statement of Theorem 1.1 is then an immediate consequence of 
(H) and (|9|). To prove the remaining estimate, let us note that it is equivalent to 
showing that if \£x\ — for every x € M, then if e > there is a A{e,p) so that 

(40) ||eA(/)||p<£A^(^')||/||2, X>A{e,p), < p < ^, 

if ex is the projection onto Vx- The endpoint p = oo follows from (12) and Lemma 3.1. 
Since it was shown in [3ol] that one has the uniform O-bounds 

(41) ||eA(/)||p<CA^(^'||/||p, P^^, 

if we apply the M. Riesz interpolation theorem (see, e.g. |So3| p. 6) to the p — oo special 
case of (|^) and (|4|) we get 



||eA(/)|l,<C"'-ei-''-A^(p)|l/|l„ A>A(e,oo), p>Hili+l), 0, = ^, 
which of course yields the remaining cases of (^). □ 

3.3. Duistermaat-Guillemin. Since it is quite simple, for the sake of completeness, let 
us point out how our pointwise bounds for the Weyl remainder term imply integrated 
ones. Specifically, if N{X) is the number of eigenvalues of \/—A that are < A, then our 
results imply that 

iV(A) = (27r)-"Vol A/A" + i?(A), R{X) = o{X''-^), 

under the assumption that |£| = 0, if £ C S*M is the set of initial directions (x,^) of 
closed loops. 
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To see this we note that our assumption |£| = and the Tonelh-Fubini theorem imply 
that one must have \Cx\ = when x G M\E with \E\ = 0. To use this, if e > is fixed, 
choose an open set fl C M containing E with < e. Since K = M\H is compact, we 
conclude from ( |T^ ) that for this s we can find a A > 1 so that |i?(A,x)| < s\"~^ when 
X G K and A > A. On other hand, we know from [ Ho 1 ] that there is a uniform constant 
C so that |i?(A,x)| < CA"-\ A > l,a; e Af. Thus, if A > A, 



\m\ 



R{X,x)dx\ < / \R{X,x)\dx- 
M Jn 



\R{\,x)\dx 



K 



< CeVol AfA"-i + eVol AfA"-\ 
which of course means that R{X) = o(A"~^), as desired. 

We should point out that the assumption \C\ = seems to be stronger than the 



assumption of Duiestermaat-Guillemin [DG| that the set of closed geodesies has measure 
zero. The two assumptions might be the same, but we do not know of a proof. 



4. Geometry of loops 



In the remainder of the paper, we apply Theorem 1.1 to a number of concrete classes 
of compact Ricmannian manifolds and investigate the extent to which it is sharp. In so 
doing, we shall need to analyze the condition \Cx\ =0 (Vx S M) and related properties 
of loops. We therefore pause in our analysis of eigenfunctions to explore the geometry of 
loops. The issues involved are motivated by analytic problems and lie someone outside 
of the mainstream of the geometry of geodesies. 

To clarify the geometry of loops, we introduce the following notions: 

Definition 4.1. We put: 

• -.^ expTiJp T*M n T*M. 

• ST^ = T'^n S*M. 

. = exp-i(x) - {TC : e e C T*M; 
. Lsp, = {|r| :r;^^0}cM+ 

• Cx e S*M : 3t e Lsp^, exp^ ^ x} C S*M 



In this list, we have mixed together notions involving the geometer's geodesic fiow (the 
Hamiltonian flow oi p^) and the microlocal analyst's geodesic flow (that of p). Thus, T'^ 
is a homogeneous closed subset oiT*M (invariant under the R+-action), and ST^ is the 
'base' of the cone. On the other hand, we deflnc exp^ in the usual geometer's, rather 
than as tTx o exptHp\T*M- For instance, we have: TSV'^ ~ {T^ : |^| — l,expxT£^ — x}. 

It is natural to say that Lspx is the 'loop-length spectrum' at the point x, i.e. the set of 
lengths of geodesic loops. The set Cx, already introduced, is the set of loop directions at 
X. The (non- homogeneous) set Tx is the set of loop (co)-vectors. In general, the sets F^, 
Tx need not be smooth manifolds. Later in this section, we shall consider restrictions on 
(Af, g) which imply smoothness properties of loopsets, and study a variety of examples. 

First, we present a standard result which plays a crucial role in several applications. 
Proposition 4.2. L* is constant on smooth submanifolds ofTx- 
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Proof. It is sufficient to prove the result for smooth curves in T^- This can be seen by using 
either the first variation formula for the lengths of a one-paxameter family of geodesies, 
or by using a symplectic argument. We present both for the reader's convenience. 

First Variation proof: Let a(s) denote a smooth curve in T^, and let 7s (i) = 
exp^tL*{x,a{s))-^^. By definition, 7^(0) = a;,7(l) = x and L*{x,a{s)) = /g Ws{t)\dt. 
Since 7s is a geodesic for each s, we have 

is lo \m\dt = Jo (£7^ m)dt 

Here, Ys{t) = -j^Jsit) is the Jacobi field associated to the variation, and ^ is covariant 

differentiation along 7^. In integrating |j by parts, the boundary terms vanished because 
^,(0) = Ys{l) = as the cndpoints of the variation were fixed. □ 

Symplectic proof: The curve C parameterized by is an isotropic submanifold of 
T*M, since it lies in T^M. Hence, the map 

i:C^ T*M, i{t^) = exp{tHp{m, 

is an isotropic embedding. Now let a = ■ dx denote the action form of T*M. Then the 
length of the geodesic loop 7s {t) (as above) is given by 

L*{x,a{s))= / ^-dx. 

J Is 

Let S denote the surface parameterized hy a{s,t) = js{t). Since all the geodesic rays 
are loops, the boundary of S equals 7si — 7so • Since S is isotropic we have, by Stokes' 
theorem, 

(42) L*{x,a{si)) - L*{x,a{so)) = f a- f ^■dx = j u) = Q. 

Here, ui is the canonical symplectic form dx Ad^. □ 

Corollary 4.3. Suppose that Tx is a manifold. If \£x\ > 0, then there exists T > such 
that TS*M c r^. 

Proof. We have just proved that the connected components of must lie in the spheres 
TS*M where £ € Lsp^. It also follows that the loop-length spectrum Lsp^ is discrete, 
hence countable. Now, is just the radial projection of T^- If it has positive measure, 
then one of the components must project to a set of positive measure, hence V'^ must 
have positive measure in TS*M. Since it is also a submanifold of TS*M, it must equal 
it. □ 

We now introduce hypotheses on (M, g) under which loopsets have good smoothness 
properties. We recall that the intersection X CiY oi two manifolds is clean if X fl y is a 
manifold and if r^(X nV) = T^X n T^Y for all x € X nY. 
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Definition 4.4. We shall say that a compact Riemannian manifold (M, g) has: 



(i) Clean loops of length T ^ ifT^ — expTHpT^M D T*M is a clean intersection. 



(a) Clean loopsets at x ifT^ is clean for all T. 

Cleanliness is related to the question of whether a one parameter family of geodesies 
at a self-conjugate point a; is a family of loops at x or not. Let us note that a curve in 
is the set of initial vectors to a one-parameter family of geodesic loops 7u(i) it G [0,T]). 
The variation vector field Y{t) — -§^^^(ftu{^) of such a family is then a Jacobi field 
satisfying Y{Q) = Y{T) = 0. Thus, each tangent vector V = ^7^1 (0) to V'^ has the form 
V = Y'ifi) where F is a Jacobi field satisfying Y{Q) = Y{T) — 0. On the other hand, any 
Jacobi field Y along a geodesic arc exp^, t£_,t < T satisfying Y{0) = Y(T) = gives rise 
to a one-parameter family of geodesies '-fu{t),t G [0, T] along which x is self-conjugate. 
However, in general 7u(i),t € [0,T] fails to be a loop at x and is just a loop to 'first 
order'. Under the cleanliness condition, this cannot happen. 

Proposition 4.5. cxpTHpT*M H T*M is a clean intersection if and only if is a 
manifold and Tf^^^^-^V'^ is isomorphic to the space Jt of Jacohi fields Y{t) along the 
geodesic arc exp^ tS,,t <T satisfying Y{Q) — Y{T) — 0. 

Proof. We obviously have 

T(x,o exp THpT*M = (d,, exp THp)Ti^.j;^^)T*M, (exp THp{x, 77) = [x, £,). 

As is well known, (d^ exp THp) is defined by 

(d„expTi7p)(J(0), J'(0)) - (J(r), J'(T)) 

where J{t) is a Jacobi field with initial conditions (J(0), J'(0)) along exptHp{x,£,),t e 
[0,T]. Clearly, expTHpT*Mr]T(^^^^)T*M is the space of value {J{T),J'{T)) of endpoint 
values of Jacobi fields satisfying J(b) = 0. It follows that r(^_j) expTHpT*MnT(^x,^)T*M 
consists of the Jacobi fields along the loop exptHp{x,£^),t G [0,T] with J(0) = = J{T). 
The cleanliness assumption is therefore Ti^xX)^x consists of all such Jacobi fields. □ 

We may reformulate the condition in terms of the length functional 

L^-.n^^M., n^{a: [0, 1] M, a{Q) = a(l) = x, a e H^[0, T]} 

on the iJ^-loop space at x. As is well-known, the critical points of are just the loops 
at X, and the set of T for which F!^ 7^ are the critical values of Lx. 

Proposition 4.6. {M,g) has clean loopsets at x if and only if Lx is a Bott-Morse func- 
tion on ^Ix i.e. a function whose critical point sets are non- degenerate critical manifolds. 



Proof. As noted in Proposition 4.2, Lx is constant on each component of its critical point 
set, so the critical point set of critical value T is precisely FiJ". To prove non-degeneracy, 
we must show that the Hessian of Lx is non-degenerate in the normal bundle of F!^ in 
T*M. The Hessian at a critical point xi is of course the Jacobi operator + i?(T, •)T 
along the geodesic through (x,^), restricted to tangent vectors to Q,x, i.e. to vector 
fields on the geodesic which vanish at both endpoints. Non-degeneracy of the Hessian is 
equivalent to the statement that its kernel lie in the tangent space to F^ , which we have 
just verified. □ 
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In the general case (without cleanUness assumptions), it follows from Sard's theorem 
(and a standard finite dimensional approximation to Q,x (Mj) that Lspx is always a set 
of measure zero in M. Putting together the previous propositions, we obtain a stronger 
statement in the clean case. 

Corollary 4.7. Suppose that {M,g) has clean loopsets. Then Lspx is a discrete subset 
o/R and is a submanifold ofT*M. 

Proof. Since = UT>or^ the second statement follows from the first and the definition 
of cleanliness. The first statement follows from the fact that Bott-Morse functions have 
isolated critical values. □ 



The following Proposition will play a role in the proof of Theorem 1.5 



Proposition 4.8. Let (M,g) be a compact Riemannian manifold. Then: 

(i) Suppose that for some T > 0, the loopset ST^ is clean and \SV^\ > 0. Then S*M = 

(a) If {M,g) has clean loopsets at x, and if \Cx\ > 0, then there exists T > such that 

Proof, (i) The hypothesis implies that ST^ is a submanifold of S*M of full dimension. 

(ii) Since Lspx is discrete, F^; = UTyoTT'^ . Since Cx is the countable union of their radial 
projections to S*M, and has positive measure, at least one component TT^ must be of 



dimension n — 1. Since TT'^ C TS*M by Proposition 4.2, again it must equal it. □ 



4.1. Examples. We now present a few simple examples to illustrate the loopset notions, 
mainly with a view of clarifying the notation and terminology. It would take us too far 
afield to provide even a sketchy discussion of the wide variety of possible phenomena that 
occur; we limit ourselves to some rather standard cases. 

4.1.1. Flat tori. Let R"/r be a flat torus. A geodesic loop at a; is a helix which returns 
to X. Loops on flat tori are always closed geodesies and they correspond to lattice points 
in F. Thus, 

F, =F, = {-^:7eF}, Lsp, = {I7I, 7 e F}, }. 

l7l 

We observe that F^, is a countable discrete subset of r^M"/F, that Cx is a countable 
dense subset of S*M, and that Tg is a countable set of embedded Lagrangean components 
each diffeomorphic to ]R"/F. Also, Ca = U^gr{(l7l, '''j a;, 0) : x G M}. This example is 
non-generic because (among other things) every loop is a closed geodesic. Flat tori are 
examples of metrics without conjugate points. We now generalize the discussion. 

4.1.2. Manifolds without conjugate points. Let {M,g) denote a manifold without conju- 
gate points, i.e. a Riemannian manifold such that each exponential map exp^ is is a 
covering map. Manifolds of non-positive curvature are examples, so this class of metrics 
is open on any M (thought it may be empty for some M). By deflnition, there are 
no Jacobi fields along any loop satisfying 1^(0) = ^(1) = 0, so the Jacobi operator is 
non-degenerate. Equivalently, each loop functional Lx is a Morse function on Qx- Thus, 
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for all countable discrete set of points, Cx is countable and is a finite set. 

Unlike the case of flat tori, loops are not generally closed geodesies. 



4.1.3. Surfaces of revolution. Surfaces of revolution provide a simple (but non-trivial) 
class of examples for which the loopsets may be explictly determined. Topologically, the 
surfaces must be either spheres or tori. It would take us too far afield to discuss the 
geometry of loops and eigenfunctions on general surfaces of revolution, so we only give 
a few indications of how to determine the loops explicitly and refer the reader to the 
articles fC^ @ pCMS| Q |z| for further background. 



A 2-sphere of revolution is a Riemannian sphere (S"^, g) with an action of by isome- 
tries. We may write g = dr^ + a{r)^d9^ in geodesic polar coordinates at one of the poles 
(fixed points) of the action. The poles are always self-conjugate points, and are points 
at which eigenfunctions attain the maximal bounds. 

Tori of revolution generalize flat tori but non-flat cases must have conjugate points 
(Hopf 's theorem) . The metric on a torus of revolution may be written in standard angle 
coordiant es a s g — dx^ + a(x)'^dd'^. Tori of revolution have no poles, and as we will see in 



Theorem 1.6 eigenfunctions need not attain the maximal bounds. In fact, although we 



do not prove it here, tori of revolution never have maximal eigenfunction growth. 

The key feature of surfaces of revolution is the integrability of the geodesic flow, which 
follows from the existence of the Clairaut integral I{x,£,) — S,{-§g)- Here, ^ is the 
vector field generating the rotations. The Clairaut integral is constant along geodesies, 
or otherwise put, unit speed geodesies must lie on level sets of / within S*M. Compact 
regular level sets are necessarily tori and are known as the 'invariant tori.' We denote by 
Tjg C S*M a. torus with Clairaut integral I = Iq- For notational simplicity, we ignore 
the fact that the level sets of / may have several components. Near any regular torus 
Ti, one can always find local 'act ion- angle' such that the geodesies may be expressed as 
winding lines on the invariant tori + tuji,(j)2 + tuj2). Here, loj are the frequencies of 
motion (which generally vary with /). 

Now consider the geodesic loops at a point x E S"^. If a; is a pole (a fixed point of 
the rotation), then all the geodesies passing through x are meridians, and the projection 
from the torus of meridians {/ = 0} to under the standard projection n : T*S'^ — > S'^ 
collapses the initial tangent vectors at x to the meridians to the point x. Thus, F^'^ — 
T*M, Cx — S*S'^, and Lspx = {nL,n G Z} where L is the length of a closed meridian 
geodesic. All of these loops are closed geodesies. 

Now let AI ~ S'^ or M — and assume x is not a pole, as must occur in the case of 
tori of revolution. Under tt : T*M M, the tori Tj project to annuli on M, and with a 
finite number of exceptions, the projection map tt/ : T/ — > M, restricted to a torus, is a 
2-1 cover with fold singularities (the caustic sets). It follows there are either zero, one or 
two geodesies of Tj passing through x, accordinly as x does not lie in ttj{Ti), lies on its 
boundary, or lies in its interior. 

Using action-angle variables, we can determine the tori Tj which contain a geodesic 
loop X. Such a loop is a winding line on T/ which passes through the set nj^lx) at least 
twice. If it passes through the same point of ttJ^(x) twice, then it is a closed geodesic; 
otherwise it is a non-smooth loop. The set of tori which contain periodic geodesies 
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are known as periodic tori, and are independent of the point x. Under a generic twist 
condition (see |Z2]), the set of periodic tori in S*M is countable. It foUows that for 



generic surfaces of revohition, the set of closed geodesies through any point is countable. 
More complicated are loops which are not closed geodesies, i.e. winding lines which pass 
through distinct points of ■Kj^ix). These do depend on x. One can write down explicit 
equations for the values of / such that T/ contains a winding line which loops at a given 
point X. It is clear from the equations that, for fixed non-polar x, a generic surface of 
revolution has only countably many loops at x. However, the question whether is 
countable at all non-polar x for generic tori of revolution app ears t o be open. For further 



information on geodesies on tori of revolution we refer to | KMS |. We plan to discuss 
this example in more detail in the future (both from the geometric and analytic point of 
view) . 

2 ^2 ^2 

4.1.4. Tri-axial ellipsoids. We recall that -Eai,a2,a3 — {ixi,X2,X3) G M.'^ '■ ^ + + — 
1} with < ai < a2 < as. Jacobi proved that the geodesic flow of i?ai ,02,03 (for any 
(ai,a2,a3)) is completely integrable in 1838. The two integrals of the motion are the 
length H{x,^) — and the so-called Joachimsthal integral J. More recent discussions 
of the geodesies of the ellipsoid can be found in 0, @, and in lCVV| (§3). 



There are four distinguished umbilic points ±P, ±Q which occur on the middle closed 
geodesic {x2 = 0}. All geodesies leaving P arrive at — P at the same time, then leave —P 
and return to P at the same time (see Q, Theorem 3.5.16 or [^). Thus, the tri-axial 
ellipsoid is an example of a Y^™-metric which is not a Zoll metric. At all other points 
X € Ea1.a2.a31 the set of initial directions of geodesies which return to x is countable p^. 

The moment map {H, J) : T* £^01,02,03 is regular away the middle closed geodesic, 

and its only singular level in S*Eai^a2,a3 is that {H = 1, J = 0) of {x2 = 0}. For a 
description of the singular level we refer to CVV[ . The middle geodesic is homoclinic, 



i.e. the trajectories on the level are forward and backward asymptotic to it. Other level 
sets of the moment map are regular Lagrangean tori, and the discussion of periodic orbits 
and loops is similar to the case of surfaces of revolution. 

4.1.5. Zoll surfaces. Now let us consider the extreme case of Zoll metrics on S'^, i.e met- 



rics all of whose geodesic are closed [Bessej. Among such metrics, there is an infinite 



dimensional family of surfaces of revolution. There is an even larger class with no isome- 
trics. 

We may suppose with no loss of generality that the least common period of the 
geodesies equals 27r. Then F^'" = 5*5^ for every x, and F^; = U^^i2TrnS*S'^ . As a; 
varies we obviously get F = U^^i2TrnS*S^. We also note that Ca is parameterized by 

R X T*M ^ T*(M X M x M), (f,x,0 ^ {t,-\£.\x,x,^,exp-tHp{x,^)). 

Since all geodesies are periodic, we obtain 

Ca = {{t, T,x,Q)} = ZxRx M, 

clearly a homogeneous Lagrangean submanifold of T*{M. x M)\0. 

As will be discussed in although geodesies are recurrent at every x, we do not 
expect eigenfunction blow-up to occur everywhere, or even anywhere in general. 
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5. Real analytic Riemannian manifolds 

Our first application is to characterize Riemannian manifolds in the analytic setting 
with maximal eigenfunction growth. 

Theorem 5.1. Suppose that {M,g) is a compact real analytic Riemannian manifold, 
with -normalized eigenfunctions {ipx} satisfying \\(/)\\\oo = ^^(A"""'^). Then {M,g) is a 
-manifold for some l,m. In particular, M = S"^ (topologically) if n — 2. 



We recall (cf. |Bessc | , Definition 7.8 (p. 182)) that a compact manifold M with a 



distinguished point m is called a 17"-manifold if all geodesies issuing from m come back 
to m at time £. We further recall that a compact 17"-manifold has the properties that 
7ri(M) is finite and H*{M,Q) has exactly one generator. 



Proof of Theorem 5.1 By assumption, there exists a constant C > and a subsequence 

of i -normalized eigenfunctions } such that | \4>\^ \ |oo > CX^ ^ . This contradicts the 
last statement of Theorem (2.1), hence there exists a point m such that //,„(£,„) > 0. 

Since g is real analytic, exp : T^M ^ M is a real analytic map, hence r„i is an analytic 
set. In any local coordinate patch U G M containing m, Tm is the zero set of a pair 
of real-valued real-analytic functions, i.e. has the form (/i(^), . . . , fniO) — • • ■ i i^n)- 
The solutions are the same as for (/i(C) — rni)^ + • • • + {fn{£.) — — 0, so is the 

zero set of a real analytic function. 

It is well-known (see [ BM , 0, y, ^ ) that the zero set of a real analytic function is locally 
a finite union of embedded real analytic submanifolds Y^"^ of dimensions 1 < fei < n — 1 . 
Thus, for each ^ there exists a ball Bs{(,) such that 

(43) T^nBsiO^^tiY'^- 

We claim that for some (^, S) there exists a component Fj""^ of dimension n — 1. If not, 
r„i is of Haussdorf dimension < n — 2. But then its radial projection p : T„i — > S^M 
would also have Haussdorf dimension < n — 2. In fact, each ray through the origin in 
T^M intersects in at most countably many points. So the radial projection preserves 
the dimension. But this contradicts the fact that = p{Tm) has positive Lebesgue 
measure. 

Now let S C r„j be an open embedded real analytic hypersurface of T*M. Consider 
the rays {t^ : < t < 1} C T^M and the union 

C = U^esK ■0<t<l}. 
Thus, each ray in C exponentiates to a geodesic loop which returns at t ^ 1. As proved 



in Proposition 4.2, it follows by the first variation formula that the length |^| of each loop 
must be a constant independent of ^ G 5. 

We conclude that |^| = for some i S K"*" and ^ €Y. But this equation is real analytic, 
hence must hold on all of S; hence S C £S^M. Again, by real analyticity, iS*^M C r„i. 
This is the same as saying that (M, g) is a Y^'"-manifold. 

□ 



In §1^, we shall prove that the converse to Theorem 1.2 holds automatically in the 
analytic setting. 
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6. Generic metrics: Proof of Theorem 1.4 



We now prove that there exists a residual subset of the space Q of C°° metrics with 
the Whitney C°° topology for which 
to prove: 



o(A 2 ). In view of Theorem 1.1, it suffices 



Lemma 6.1. There exists a residual set TZ G Q of metrics such that \C^\ = for every 
X ^ M when g £ TZ. 



Proof. Recall that \C^\ = if and only if /g^^^ l/L*{x,S,)d£_ = 0, if is the loop- 

length function defined before. 

To make use of this, choose a coordinate patch fl C M with coordinates y = k{x) 
ranging over an open subset of R". We then K C K(f2) be compact and let 



F{g) = sup / 
yeK Js" 



da 



using the induced coordinates {j/,^} for T*r2 C T*M. Here also, da is the standard 
surface measure on 5""^, and we are abusing notation a bit by letting L*{y,^) denote 
the pushforward of L* using k. It then suffices to show that the set of metrics for which 

GN^lg: Fig) < 1/N} 

are open and dense. 

Density just follows from the fact that F{g) — for any non-ZoU real analytic metric. 
Such metrics are dense in Q. 

The main step in proving that these sets are also open is to show that the function 
f{9iy) = ^sn-i dcf I L*g{]j,i) is upper semicontinuous on X ii'. This holds since 1/L*(y, ^) 
is a positive, (locally) bounded upper semicontinuous function on ^ x x 5"^^, if we 
equip Q with the C^-topology. Therefore, if (gj^yj) —> {g,y), Vj G K, we have 

f da f da 



, f da f da f da 

using the dominated convergence theorem and upper semicontinuity. 

Now let us prove that the sets ^at are open. Let g G Qn. By definition of F, f{g, y) < 
for each y £ K. Since / is upper-semicontinuous, the set {/ < -^1 is open, so there exist 
6{y) such that Bg(^y-^{y) x B^i^y-^{g) C {/ < ;^}, if Bs{y) and Bs{g) denote the 5 Euclidean 
and balls of y and g, respectively. Here, g is fixed so we do not indicate the dependence 
of the (5's on it. As y varies over the balls Bs(^y-^{y) give an open cover of K , and by 
compactness there exists a finite subcover {Bs(^y.-){yj), j — 1, . . . , N}. Let 6 = miuj S{yj), 
so that Bsiyj) x Bsig) C {/ < ^} for j = 1, . , N. If .9' e Bs{g), then fiy,g') < i for 
all y G K. Hence F{g') < and so Qn is open. □ 
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7. Converse results: Proofs of Theorems 1.5 and 1.6 



The purpose of this section is to prove Theorems 1.5 - 1.6. We shall need to go into 
somewhat more detail about loopsets. 

We first consider the simplest converse in which the loopsets are assumed to be clean. 



7.1. Theorem 1.5 for manifolds with clean loopsets. 

Theorem 7.1. Suppose that {M,g) is a compact Riemannian manifold. Then: 

(i) //r^ ^% is clean, and {T} is isolated in Lspx, then R{X,x) — o(A"^^) =4> |5T^| — 
0. 

(ii) // (M, g) has clean loopsets. Then 

(44) i?(A,a;) =o(A"-i) =^ =0. 

We begin the proof with some considerations which are valid for all compact Riemann- 
ian manifolds. 

For each (T, x) we introduce the sequence of probability measures 
(45) 



A„<A 



on S^. Here, S^iTx^ (6) denotes the point mass at the indicated point of 5*^. 

Lemma 7.2. Let {M,g) be any compact Riemannian manifold. Suppose that R{X,x) 
o(A"^^). Then, ht,x,\ dO for all T . 



Proof. The fcth Fourier coefficient of fJ.T,\.x{S) equals 
(46) 



Hence R{X,x) — o(A" ^) implies 



(47) 



Ex{x,x)fiT,x,4S) - Jo e'*^^dEx{x,x) 

= e'^^^dX" + e'^'^^dR{X, x) 

= nj^ e''=^^A"-MA + e'^^^i?(A, x) - ikT e^'^^RiX, x) 



= o(A"), for k^O 

Since d9 is the unique probability measure /i on S^ satisfying fl{k) = Ska, we have 



(48) 



ike 



SI 



dfiT^xAO) ^0 (Vfc ^ 0) ^ {dfiT.x,x{0) ^ dd}. 



Hence Lemma 7.2 follows from (|47|)-(48). 



□ 



In view of Lemma 7.2, to complete the proof it sufhces to show: 



26 



CHRISTOPHER D. SOGGE AND STEVE ZELDITGH 



Lemma 7.3. We have: 

(i) Suppose that {T} is an isolated point in Lspx and that F!^ is clean. Then ^it,\.x{&) ^ 
dO implies that SV"^ is a submanifold of S*M of dimension < n — 2. 

(ii) (itaAO) ^ de (VT^O) =^ \Cx\=0. 

Proof, (i) By (^), (i) asserts that if is clean, then 

(49) / e''''>dtiT,xAd)^0 (Vfc) =^ |5r;^|=0. 

We prove this by studying the asymptotics of [U{T)E\]{x,x) as A — > oo, which are 
dual to the singularity of U{T + s){x,x) at s = 0. These asymptotics were studied (in 
greater generality) in j^] ( §2), and wc follow that discussion in the calculation to follow. 

Following [||, we write t/(r+s,a:,a:) = J^^^^ U{T+s,y,y')6x{y)Sa;{y')dydy' . Viewing 
X and T as fixed, we obtain a distribution in s with singularities at times in the loop- 
length spectrum Lsp^ (referred to as 'sojourn times' in the more general context of |zj). If 
{T} is isolated in Lsp^, then there exists an interval around s = so that U{T + s, x, x) 
is singular only at s = 0. If T'^ is a clean intersection, it then follows from the clean 
composition calculus that U{T + s, x, x) is a Lagrangean distribution in s in this interval 



(see §2, for the specific case at hand, or |DG] [Ho IV] for the general theory). This 
implies that U(T+s, x, x) has an isolated singularity at s = with a singularity expansion, 

U(T + s, x, x) = OriT, x){s + iOy + a,._i(T, x){s + iO)-''+'^ + ■■■ 

where • • • denotes smoother terms, and where the order r — dimSTj. Equivalently, we 
may write 

(50) U{T + s){x, x) = r ^M^e-^dA, 



where At{x, A) is a symbol of order r — 1 in A. By Lemma 2.5 or by (j^]. Lemma 3.1), 
it follows that 

(51) [U{T)E):\{x, x) - At{x, A) = 0{\^~^). 
Therefore, the statement (i) in Lemma |7.3| is equivalent to 

(52) order of At{x,\) <n-l ^ \TI\ = 0. 

The principal symbol of U{T + s, x,x) a,t s = (for fixed x) was calculated in (jzj. Propo- 
sition (1.10)-Lemma(3.6)-Theorem (3.15)). The formula depends on dimST^J". When 
dim5r!^ — n — 1, then the order of At equals n and the symbol is given by 

(53) aniATix,-))^\Srl\. 

Now suppose that dfiT.x,\ ~* d9. Then the first moment of dfiT,x,x tends to zero. If 
dimST^ = n — 1, it must equal S*M. But then the order n symbol of At cannot equal 
zero. This contradiction completes the proof of (i). □ 

(ii) This ca se f ollows from (i) and is easier since the hypothesis is much stronger. By 



Proposition IIJ, if \Cx\ > 0, then there exists T > such that = r^*^^. As just 
proved, this implies that the symbol of order n of At is positive. But R{x, A) = o(A~2— ) 
implies that it equals zero. □ 
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7.2. Real analytic metrics. We also have a converse theorem for real analytic metrics, 
which is a simple corollary of the above theorem. 

Corollary 7.4. Let {M,g) be real analytic. Then R{x,\) = o(A""i) =^ = 0. 

Proof. The set of vectors e S** : exp^ ^ = x} is an analytic set. Hence, it has 
only countably many components. Since L* is constant on smooth components, L.spx is 
countable for each x. 

Suppose now that \Cx\ > for some x. Then, (as discussed in the proof of Theorem 
|, there must exist T > such that ST^ = TS*M. The time {T} is isolated in Lsp^. 
Moreover, is necessarily a clean intersection, since by homogeneity oiexpTHp we 
have {expTHp)T*M = T*M. Hence, the estimate follows from Theorem |t[ □ 

7.3. Almost Clean. We now generalize the result to 'almost-clean' loopsets. In the 
argument above, we assumed that for (T, x), U{T + s){x, x) is a Lagrangean distribution 
in s, but much less is sufficient since we only used the principal symbol and we only 
required o(A") from the remainder estimate. 

Definition 7.5. We say that T'^ is an 'almost- clean' loopset if, C T*M is a (ho- 
mogeneous) suhmanifold with boundary and if the intersection expTHpT*M nT*M is a 
clean intersection in the interior of . We say that {M,g) has almost clean loopsets if 
this property holds for all T. 

Simple examples of metrics with an almost clean loop-set are provided by generic 
bumps on standard spheres. In other words, let Bx{r) denote the ball of radius r around 
a point X € (S"", go) (the standard sphere), and let e^"(?o denote a conformal perturbation 
with supp u = Bx{r). For points x in the exterior of the 'bump' B^ir), V"^ contains any 
direction which avoids the bump. For generic u it will not contain any other loops at 
time T — 2tt. Hence, S'F^'^ will be a submanifold with boundary of S*M, with boundary 
the set of initial tangent vectors to geodesies which intersect dB^ (r) tangentially. More 
generally, we could put several bumps on the sphere and obtain a number of components. 
Instead of bumps, one could put handles, and thus obtain any topology. This example 



will be important in the proof of Theorem 1.6 



Since could have several components, there is a possibility that cancellation between 
components can cause slow growth in U{T)E\{x,x) even when \£x\ > 0. We adopt a 
reasonable hypothesis, easy to check in our main examples, which rules this out. 

Theorem 7.6. (i) Suppose that {T} is an isolated point in Lspx of {M,g), that F^ is 
an almost-clean loopset , and that the (common) Morse indices of loops in all components 
ofT^ of dimension n — 1 are the same. Then R{X,x) = o(A"~^) =J> |»S'F!^| — 0. 

(ii) Suppose that {M,g) has almost-clean loopsets, and that, for all T, the (common) 
Morse indices of loops in all components of of dimension n — 1 are the same. Then 

(54) i?(A,a:) -o(A"-i) =^ \Cx\=0. 

Proof. We only prove (i) since, as in the clean case, (ii) is an immediate consequence. 



Nothing changes from the proof in the clean case until (|50|). Since the clean compo- 
sition hypothesis is not assumed to hold, U{T -\- s,x,x) is not known to be Lagrangean, 
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although the singularity at s = is still assumed to be isolated. We therefore cut the 
wave kernel up into pieces which can be controlled easily. The following argument is 
somewhat similar to the proof of Theorem (3.20). 

We fix {T,x) and choose a cover S*M adapted to ST^ as follows. By assumption, 
STJ is a submanifold with boundary of S*M for all T. It therefore has only finitely 
many connected components. We write 

where {K^^} are the components (possibly empty) of dimension n — 1 and {N]^^} are 
the components of dimensions < n—2. The Morse indices of the loops in each component 
if J-' are the same; we are assuming the Morse index is the same constant for all if J-*, as 
is the case in many natural examples. 

We then form the cover {intK]^^ ,T^{dK^^),T,{N^^), O}, where T,{dK^^) is the tube 
(tubular neighborhood) of radius e around the boundary of if J-', where T^{N'^^) is the 
e-tube around iVj^ , and where O covers the rest of S*M. 

We then define a partition of unity subordinate to the cover as follows: we define 
functions € {intK^^) which are invariant under expTHp and = 1 in points of 
miifP' which he outside of T^{dK^^). We further define 4 G C§°{T^{dK^^)), G 
C^{T,{NP)), and r, € Cg°{0) so that EJIi K + EJIi 4 + E^Ii di +r, = l. We note 
that 6, c, d, r also depend on T; since it is fixed throughout the calculation, we do not 
indicate this in the notation. 



In a well-known way (see [ Ho IV j), we may quantize the symbols to obtain a pseudo- 
differential partition of unity: 

ut tlt Tt 

j=i i=i i=i 

We then have 



UiT)E^{x,x) = [J2BlUiT)Ex]{x,x) + [J2CiU{T)Ex]{x, x) + [Y^DlU{T)Ex]{x,: 
j=i j=i j=i 

+ [R,U{T)Ex]{x,x). 



Since is microsupported in the set where V'^ is a clean intersection, it follows that 
BIU{T -\- s)\{x^x) is a Lagrangean distribution for s near and its principal symbol of 
order n is given by (i™^ /gr^ bidvx) as in the clean case. Here, rrij is the Morse index 
of loops in if J-' (which we a ssume independent of j). It follows by the standard trace 
asymptotics (see Lemma ^.5[ ) that 

[BlU{T)Ex]{x, x) = z'"^- ( / bldi^x)X" + 0(A"-i). 
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Also, R^U{T + s,x,x) is smooth is smooth in s g M since no loops at x occur in the 
microsupport of R^. We conclude that 

(55) U{T)E).{x,x) = (V / bldy^)X' + ([V CiU{T)Ex]{x,x)) 



+ [Y,DlU{T)Ex\{x,x)) + 0{X-'). 

We estimate the two middle terms in the same way; the harder one is the C-term, so 
we leave the details of the D - terms to the reader. We have: 

\CiU{T)Ex]{x,x)\ ^\T.K<xe^^-^'M^)CiM^)\ 



(56) 

= 0{^Vol{T,{dKl')X') = 0(V^A"). 
In the above we used that, for any 0th order pseudodifferential operator, 

V |C,0.(x)|2 = CExC*{x,x) - ( / |c|2dz.,)A" + 0(A"-i). 

Since T^{dK^^) is a tube around a hypersurface in 5*Af, the integral is 0(e). Thus, by 
the same Fourier coefficient calculation as in the clean case, R{X,x) — o(A"~^) implies 
that U{T)Ex{x, x) = o(A") for aU T 7^ 0. It follows for any e > that 



By the assumption that rrij = m for all j, it follows then that jST^j = 0. This completes 
the proof of (i). 

To prove (ii), it suffices to observe that Lspx is countable, and we have just seen that 
each Sr'^ has measure zero. □ 



7.4. Counterexamples: Proof of Theorem 1.6. We now present an example of a 
C°° (but not analytic) surface (M, g) satisfying: 

3a: : |£^| > 0, |i?(A,a;)| = f7(AT^) but L°°{\g) = o(At^). 

Thus, neither \Cx\ > nor |i?(A,x)| — is sufficient to imply maximal eigenfunc- 

tion growth. 

Our example is of a surface of revolution, i.e. a surface (M, g) whose metric is invariant 
under an action of by isometrics. The two possible M are M = S'^ and M = T"^ (the 
two-torus). On AI — we can construct metrics such that 3a; : \Cx\ > 0, but no 
sequence of eigenfunctions has maximal growth at x; but the invariant eigenfunctions 
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will have maximal growth at the poles. To obtain an example of a surface which does 
not have cigcnfimctions of maximal growth, we instead choose M = T^. 

We often write M = x with the angle coordinates {x,9). We shall refer the 
circle over which x varies as the base circle, and the orbits as the fiber circle. Our torus 
of revolution is of the following kind: 

Definition 7.7. (T^, g) will be called a torus of revolution with an equatorial band if in 
the standard coordinates {x,y) on we have Qa = dx^ + a{x)^d6^, where a = Vl — x'^ 
for X G (— £, e). 

Such a surface may be obtained by by revolving the graph y = e^^^\ < a; < 27r of a 
27r-periodic strictly positive function e"^^^ around the x-axis, and then joining the ends 
X = —l,x = 1 together. The 'equatorial band' is the part ~ {^i^) ^ lying over 
angular interval (—£,£) (for some < e < 1/2), and is isometric to an equatorial band 
lying over the same interval of the ^;-axis of the standard S"^. Clearly, |£^| = 27r for the 
open annulus of x G B^. 

We now construct u such that the eigenfunctions of {T'^,gu) fail to have maximal 
growth. We break up the base circle x G into four parts: a flat part on [—1, — 2£] U 
[2£, 1], a spherical part [— £, e] and a 'bridge' over (— 2£, — £) U (£, 2£) between them. We 
shall construct a metric which is invariant under the involution x —x (i.e. e'^ e~'^), 
as follows: Introduce C°° cutoff functions on satisfying: 

{1, XG [-£,£] 
0, a;e [-1,-2£]U[2£,1] 

{0, X G [— 
1, [-1,-2£]U[2£, 1] 

We now define a class of profile curves of the form: 

af{x) = a/1 -a;2xi(a;) + - s'^X2{x) + f{x){l - xi{x) - X2{x)), 

fix) = f{-x), /GC°°([£,2£]). 

For concrctcncss, we shall constnict / so that a/ is monotonically decreasing from 
the boundary of the round region (where a/(£) = Vl — e^) to the fiat region (where 
a/ = — £2). We shall begin with one such function /o and then consider small per- 
turbations f — fo + £u. We denote the associated surface of revolution by {T'^,gf), and 
its Laplacian by A/. It is clear that only the values of / on ±([£, 2£]) will matter. We 
equip the space of m G C°°{[e, 2e]) with the Whitney C°° -topology. 

Our first result involves the splitting of eigenspaces Ex = ffi^'^'f"^' Va,„ of A/q into irre- 
ducibles for S^. Here, V\^n denotes the joint eigenspace of A, Because the eigenvalues 
of Sturm-Liouville operators on have multiplicity at most two, and dim Va.h = 2, the 
multiplicity of each Vx in Ex is at most one. Therefore, mult{X) counts the total number 
of distinct irreducibles in Ex . The following Lemma shows that by a generic small pertur- 
bation of /o, we can arrange that all of the eigenspaces are irreducible. We note that the 
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same proof shows more generally that all eigenspaces are irreducible for generic surfaces 



of revolution. This result is analogous to the earlier result of one of the authors [Z3| that 
all eigenspaces are irreducible for generic G-invariant metrics on Riemannian d-manifolds 
when G is a finite group and d > do for some minimal dimension dc depending on G. 



We refer to |Z3] for background, references and arguments that are valid when G = S^. 



Lemma 7.8. Let Mirr denote the subspace of C°° {[e , 2e]) of functions f such that all 
eigenspaces of Af are irreducible for the action. Then Adirr is residual in G°°([£, 2e]). 
In particular, its intersection with any closed ball Bs{fo) around fo in C°^'{[e,2e]) is 
residual in the ball. 



Proof. We need to prove that Adirr contains a countable intersection of open dense sets. 



As in |Z3|, we denote hy A4k the metrics for which all of the first k eigenspaces are 



irreducible. It suffices to prove that each A4k is open-dense. 

Since eigenvalues and eigenfunctions move continuously with perturbations, it is a 
simple fact based on the minimax characterization of eigenvalues that Alfe is open in 
C°°{[e, 2e]) and hence that its intersection with any closed ball is a relatively open subset 



of the ball. We refer to [Z3| for the standard argument. The novel point to prove is the 
density of each Mk- 

We now separate variables. For each n e Z, we denote the joint eigenvalues and 
normalized eigenfunctions of A, ^ by {A^ $„j}. Thus, $„j(a;,0) = e™^<i>„j (x, 0). We 
shall refer to such joint eigenfunctions as 'equivariant eigenfunctions'. With no loss of 
generality, we may conjugate A by ^/a to put it into Liouvillc normal form (i.e. transform 
it to the half-density Laplacian). The functions {x, 0) are then complex eigenfunctions 
of the associated Sturm-Liouville operator 

Ij -\- + V ^^/"^ ^ ~^ ^ 

" dx'^ a(r)2 ' dx dx 



When n 7^ 0, we have 



<i>„,\ix,e) = e"'''{unA^) + iv„,x{x)), 



where Un^x,Vn^x are the real eigenfunctions of L„. When n = 0, there is only a one- 
dimensional space of invariant eigenfunctions. 

To prove that is dense, it suffices to prove that for each of first k eigenspaces, it is 
possible to find a perturbation which splits the eigenspace into irreducibles. Since only 
a finite number of irreducibles are involved, it suffices to do so for one eigenspace at a 
time and for one pair of irreducibles in its decomposition. The case where n = for one 
of the eigenspaces is somewhat different, so we assume n, m 7^ at first. 

We argue by contradiction. Suppose that, for a metric Qq, the eigenspace Ex contains 
at least two complex normalized eigenfunctions ^n,XT^m,x with m,n 0. We would 
like to find a curve of metrics which splits up their eigenvalues. It is clear that m ^ n 
since the eigenvalues of Sturm-Liouville operators are at most double. We recall that the 
variation of the eigenvalue under a deformation e^"=f/o is given by 
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where A is the variation of A. The same calculation holds for ^m,X- If the eigenvalue 
does not split under any deformation, then we have 

{u^n,X, $n,A) = (M*m,A, *m,A), (V« G C°°(e, 2e)) 

Since ^ 

/ (|$„.A(x,0)|2-|«'„,A(x,0)n7idx = 0, (VueC°°(e,2e)) 
Jo 

it follows that \^n,\{x,9)\ = \'i/m,x{x,0)\ in (£,2£) x S^. It follows that on this region 
^n,\{x,0) = e'("~'"^^e'"(^)'J'm,A(a;,^) for some real-valued function a. 

To analyze the situation, we write 
Then Ae*"^e^(^)+*'^^^^ = Ae*"^e''(^^+*'^^'^^ gives 



/ 2 

a' , , 



p" + zr" + iip'y + 2it'p' - {r'f) + -{p' + it') = A + ^ 



a 



Again separating into real and imaginary parts this becomes: 



(57) 



[ 9: = 2p't' + t"{x) + ^t'{x). 
Prom the imaginary part equation we obtain 

(58) / = e^'W(a-ie-20 =^ r = t(0) + e^'(°) / {a-^e-^P)dy. 

Jo 

Thus, the phase r is determined up to two constants t(0),t'(0) by the modulus of 
We further observe by combining this with the real part equation that p solves: 

/ 2 
J' I I ^ ' ^2-r'(0)/„-2^-4p\ ^2 

a a[x)^ 
Suppose now that pi = |<I'|,P2 = |^|- Subtracting their two equations gives 

Here, we assumed that ti(0) = T2(0), which is no loss of generality since we may multiply 
\1/ by a e^^'^'^^^^~'^^^^^^ obtain a scalar multiple with this property. 

If we put f = pi — P2, then we obtain the inequality 

£f = ip'l + P'2 + ^).f' «-'e2-(0)/G(pi, P2) 

=^ \f"{x)\<A\f{x)\+B\f{x)\ 

for some constants A, B, namely the suprcmums of (p[ + P2 + ^) and G respectively. It 
follows then by Aronsajn's unique continuation theorem that / = on an interval implies 
/ = 0. Hence we have that \^{x)\ = \^{x)\ globally. 

Of course, this is no contradiction: there do exist metrics on possessing pairs 
equivariant eigenfunctions with the same eigenvalue and same pointwise norm: indeed. 
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for the standard E^/^^ one may take $„,A(a;,6') = e*('=^+"^), *„,A(a;, 6*) = e'(^^+™^), 
where fc^ + = + We now assert that such pairs only exist for flat metrics of 
revolution on . 

Indeed, suppose that $ = F*. We also note that F{x,e) = e^(""")''F(a;, 0). Write 
F(x, 0) = /(a:), ^(a;, 0) = i'ix). Since also |/(x)| = 1, there exist smooth functions p, r, a 
such that 

^X^F^ = A(F*) = *AF + ■ VF + FA^ 4=^ -^AF + • VF = 

^=4> V(^'2VF) = 0. 

Simple calculations give: 

V(*2VF) = =^ ^/(a:) - 2(f£^/'(x) + (fix) + 

=^ ^ = 2{p'{x) + ^T'{x))[^a'{x)] - a' {xf + ^a" {x) + 
Breaking up into real and imaginary parts gives 

^ = ^2T'{x)a'{x)-{a'{x)f 

3: Q^2p'{x)a'{x)+a"{x)-^a'{x). 
The imaginary part equation has the unique solution 

a'[x) = (n-m)a-i(x)e-2p(^) 
consistent with the given initial conditions. Substituting into the real equation gives: 

(59) m^~n^ = -2{n - my {x)a-\x)e-^P'^''''^ - {n - mfe-^P'^''\ 

Substituting ( p8| ) into ( |59| ) gives 

(60) m2 - n2 = Ce-^P => p{x) = Cn,™ 

for some constant Cm,n- Substituting this into the real part equation of ( p7[ ) and substi- 
tuting the value of r' in (|5^) gives 

(61) [C^ - r?)a-^ ^ X, 

which implies that a{x) is constant. □ 

To complete the proof we use the following Lemma. 

Lemma 7.9. For any 'bridge function' f , the joint -normalized eigenfunctions {^„ \} 
of A and of the S-^ action on {T'^,ga) satisfy ||$n,A||oo — o(A~2~). 

Proof. Let us now prove eigenfunction estimates on our torus of revolution T with metric 

(62) ga = dx'^ + {a{x)Yde'^, -1 < a; < 1, -tt < 6* < tt, 

where of course a is a positive periodic function on K/Z. We shall show that for L^- 
normalized eigenfunctions of the form 

(63) $ = $„a = e*"^$(x) 
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with eigenvalue A we have the bounds 

(64) miL-iT^) < 

for some cr > to be specified later. These bounds of course imply that for our the above 
metric we have L°°(X,g) = o(A^/^); however, as we shall sketch after the elementary 
proof, the bounds we obtain from this simple argument are far from optimal. 

To prove (^J), fix an even function p e C^(R) with integral one and having the 
property that 

p{t) = if 1^1 > 1. 
Then, if p is the Fourier transform of p, it follows that for a given i5 > 

(65) ■!> = xi* = ^(^(A - y^))$. 
Note that 

(66) = 2^ / ^-''^'--^^pis-'t) at. 

Let x^(w, z), w,z £ , be the kernel of X\- Then by finite propagation speed for 
cos{t\/ — A) it follows that 

xiiw,z) =0 if dist(it;,2:) > 26, 

with dist(it;, z) being the distance between w and z measured by the metric ga- Because 
of this, we can explain the role of S in our estimates. First of all, if (5 > is small 
the kernel is supported on small scales, and so after rescaling to the unit scale in local 
coordinates behaves increasingly like ones arising from a flat metric as S ^ 0. This 
makes calculations in the proof easier; however, as (5 ^ one sees from (^5|) that X\ 
looks less and less like an eigenspace projection operator, which will reflect negatively on 
(U). Ideally, one would like to take S to be, say, half the injectivity radius, but then the 
calculations that arise in the proof seem formidable. 

To proceed, let us note that if S is larger than A^""", where ao > is fixed and if 6 is 
smaller than half the injectivity radius then we can write 

± 

where the amplitudes satisfy the bounds 

with C-y independent of 6. This follows from routine stationary phase calculations using 



(|66|). (The case 5 « 1 was carried out in |So3|, p. 138-140.) Also, the support properties 
of Xxiw, z) imply that 



g{w, z) ~ if dist(ii', z) > 2S. 

To p roceed, let us fix w € and use the natural coordinates {x, 9) corresponding 
to (|6^. In proving bounds for (f> sd w, in view of (^3|), we may assume that w has 
coordinates (a:o,0). If {x,9) then are the coordinates of z, it follows from ( |63|) that 

(67) 0(a;o,O) = ^ J K^^^^g{xo, x)(j){x) dx, 
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where g is given by the osciUatory integral 



By applying the Schwarz inequality we of course get 



± 



using in the last step the fact that dVol a{x)dxd6 with a{x) > and the fact that 
$(a;, 0) has L^-norm one. 



Because of this our proof of boils down to estimating the norms of K^^ „ ^ 
and then optimizing in 5. To do this note that the integrand of the oscillatory integral 
defining ^ vanishes unless \9\ + \xq — x\ < C5. Therefore, it is natural to make the 
change of variables 9 — > 59 ^ a; — > .tq + which results in 

(68) ifA,„,^(^0' ^0 + = i^l^)^'^ j e±*^'*^(""^^''^)at5(^"' 0; + 
where 

= 5" Mist ((xo, 0), (a;o + 59)) + n59. 

If the warping factor a{x) for the metric (^) were constant in a 26 neighborhood 
of xq then one would have {d/d9)^^s = -{a{xQ)yf /\y'^ + {a{xQ)9f\^/'^, which results 
in a lower bound of (a(a;o)y)^ for the Hessian of the phase function in this case. For 
non-constant a one has 



5-Mist((a;o, 0), [x^ + 5y, 69)) ^ y/y^ + {a{xo)9)^ + rs{x„,y, 9), 
on the support of the oscillatory integral where for where for a given j — 0, 1, 2, 

\{d/d9)^rs\ < C6. 

Consequently, there are uniform constants cq > and Cq < oo so that 

\id/d9)Hs\>coy^ if y^ > Co5. 



Since the amplitude in (|68|) has uniformly bounded derivatives, we conclude from 
stationary phase that 

\K^^^Axo,xo + 6y)\<C{X/5)^^'{yHX)-'^'^C6-^\y\ if \y\>Co5-'^\ 

One also has the trivial bound l^^^^^l < C{\/6)^^^ . Using these two bounds we can 
compute the norm of the kernel K^^^ with respect to the original coordinates: 

(/ \K^,„Axo,x)\'dx)'/'=5'/'{l \Kl^^,{xo,xo + 6y)\'dy)'^' 

<C6-^I\\ y-'dy)'/' + CX'/^{ dy)'^' 
Js^/^ Jo 

< C6-^'^ + C\^/H^/\ 
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One optimizes this by choosing S ~ X which impUes that one has (^) with a = 
1/8. □ 

Completion of the proof of Theorem 

Two of the three properties follows from the corollary and from the fact that \Cx\ > 
within the band. For the third property, we need to prove the existence of a point xq such 



that i?(A, a:o) = fl{X^~). For this we appeal to Theorem 7.6 (i). Let a; be a point on the 
equator of the round part. Then F^'^ has an almost-clean contribution coming from the 
great circles in the equatorial band. We claim that by choosing the flat part to be long 
enough, there T'^^ contains no other loops. Indeed, any other loop would have to enter 
the flat part. It cannot reverse direction, so it must travel around the flat part before 
re-entering the spherical part. As long as the flat part is at least 27r units in length, any 
such geodesic loop at x will not belong to T^^ . □ 

We believe that the correct bound should be X^^^. We shall discuss this further below. 

8. Further problems and conjectures 
Let us mention some related natural problems which remain open. 

Problem 1 First, what is a sufficient condition for maximal eigcnfunction growth at a 
given point x or at some point of {M,g)l 

As mentioned above, > is not sufficient to imply the existence of a sequence of 
eigenfunctions blowing up at x at the maximal rate. In fact, to our knowledge, the only 
Riemannian manifolds known to exhibit maximal eigcnfunction growth are surfaces of 
revolution and compact rank one symmetric spaces. In the case of surfaces of revolution, 
invariant eigenfunctions must blow up at the poles because all other eigenfunctions vanish 
there and yet the local Weyl law must hold. In the case of compact rank one symmetric 
spaces, the exceptionally high multiplicity of eigenspaces allows for the construction of 
eigenfunctions of maximal sup norm growth. In each case, Cx has full measure, but 
the mechanisms producing maximal eigcnfunction growth involve something more. All 
of these examples are completely integrable and eigenfunctions with maximal sup-norms 
may be explicitly constructed by the WKB method. Let us describe the symplectic 
geometry, because it turns out that we are tantalizingly close to proving it must occur in 
the real analytic case. 

The eigenfunctions with maximal sup norms in these cases are actually oscillatory 
integrals (quasimodes) associated to (geodesic flow)-invariant Lagrangean submanifolds 
diffeomorphic to S'*M x S"! - S"""! x 5^ They are the images of S'^ x S';„M under the 
Lagrange immersion 

i-.S^x S*^M T*M\0, = G^{x,^), 

where G* : T*M\0 — > T*M\0 is the geodesic flow. Under the natural projection tt : 
S*M M, TT : i{S^ X S;^M) M the sphere S*M is 'blown down' S*M to m. As 



discussed in |TZ] (and elsewhere), singularities of projections of Lagrangean submanifolds 
correspond closely to sup norms of the associated quasimodes, and such blow-downs give 
the maximal growth rate of the associated quasimodes. Existence of a quasimodes of 
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high order attached to an invariant Lagrangean S'^M x may therefore be the missing 
condition. 



In the real analytic case, we show in Theorem 5.1 that all geodesies issuing from some 
point m return to to at a fixed time i. If, as is widely conjectured, they are smoothly 
closed curves at to, then these geodesies fill out a Lagrangean submanifold of the form 
b{S^ X S^M) above. To complete the conjectured picture, we would to show that there 
exist eigenfunctions which are oscillatory integrals associated to this Lagrangean. It is 
automatic that quasimodes (approximate eigenfunctions) can be constructed, but it is 
not necessarily the case that they approximate actual eigenfunctions. 

We do not expect such quasimodes always to approximate eigenfunctions, nor do we 
expect maximal eigenfunction growth in all situations where Cx has full measure, even 
at all points x. In other words, we do not expect maximal eigenfunction growth on all 
Zoll manifolds (manifolds all of whose geodesies are closed) . It is quite conceivable that 
(non-rotational) Zoll spheres do not have maximal eigenfunction growth, even though the 
converse estimate i?(A, x) — S1(A"^ ) holds. Indeed, it is known (Q) that on the standard 
S^, almost every orthonormal basis of eigenfunctions satisfies ||^a|| = 0{^/\og A), even 
though special eigenfunctions (zonal spherical harmonics) have maximal eigenfunction 
growth. It is possible that eigenfunctions of typical Zoll surfaces resemble such typical 
bases of spherical harmonics rather than the special ones with extremal eigenfunction 
growth. In short, the converse direction appears to be a difhcult open problem. 

In a future article with J. Toth, we plan to study L^-norms of eigenfunctions and 
quasimodes directly using oscillatory integral and WKB formulae. As a very special 
case, we will improve the sup norm bound on eigenfunctions of tori of revolution to A^/^. 

At the opposite extreme is the problem of characterizing compact Riemannian man- 



ifolds with minimal eigenfunction growth, such as occurs on a flat torus. In |TZ| it is 



proved that in the integrable case, the only examples are flat tori and their quotients. 
Problem: Characterize (M, g) with maximal L^'-norms of eigenfunctions. 

We have left this problem open for p < ^^^z-p-, and we expect the condition on (M, g) 

to change at the the critical Lebesgue exponent p — ^^^[^p-- Indeed, as one of the authors 
proved in [ |So2| , the 'geometry' of extremal eigenfunctions changes at this exponent. 
Specifically (cf. [3o3|, p. 142-144, [3ol|), forp > '^'"^^P eigenfunctions concentrated near 
a point tend to have extreme norms, while for 2 < p < "^^^^P ones concentrated along 
stable closed geodesies tend to have this property; if p ~ ^''^-P ' least in the case of 
the round sphere, both types give rise to f](A^(P)) bounds. 

Bourgain Q has constructed a metric on a 2-torus of revolution for which the maximal 
bound is attained although \£x\ = for all x. The eigenfunctions are similar to the 
highest weight spherical harmonics on 5*^ which concentrate on the equator. Note that 
this is the Lebesgue exponent where one expects the behavior of eigenfunctions which 
maximize this quotient to change. 

Thus, we do not expect > 0' to be a relevant mechanism in producing large 

LP norms below the critical exponent. In some sense, existence of stable elliptic closed 
geodesies is more likely to be involved. We plan to study the problem elsewhere. 
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